Theoretical model of donor–donor and donor–acceptor energy transfer on a nanosphere

In this study, we introduce a novel advancement in the field of theoretical exploration. Specifically, we investigate the transfer and trapping of electronic excitations within a two-component disordered system confined to a finite volume. The implications of our research extend to energy transfer phenomena on spherical nanoparticles, characterized by randomly distributed donors and acceptors on their surface. Utilizing the three-body Padé approximant technique, previously employed in single-component systems, we apply it to address the challenge of trapping within our system. To validate the robustness of our model, we conduct Monte Carlo simulations on a donor–acceptor system positioned on a spherical nanoparticle. In particular, very good agreement between the model and Monte Carlo simulations has been found for donor fluorescence intensity decay.

decay or emission anisotropy decay.Therefore, theoretical models traditionally applied to describe disordered infinite systems in terms of energy migration and transfer may not be universally applicable to nanostructures of finite volume, specific geometries, and limited fluorophore numbers.
In the our previous works 28,29 , we discussed energy transfer in spherical nanoparticles with a finite radius R , where molecules attached to the surface act as energy donors.In our new paper, we analyze spherical nano- particles where both donors and acceptors are attached to the surface, creating a more complex system with additional physical processes.The development of a specialized theoretical model capable of elucidating these intricate processes would not only fill a critical gap in our understanding but also pave the way for the precise design of nanoparticles with enhanced functionalities.From advanced imaging techniques to targeted drug delivery systems, the implications of such a model span a wide range of cutting-edge applications, promising breakthroughs in both research and practical implementations.
In this investigation, we introduce a novel approach utilizing the three-body Padé approximant method applied to the Green function, developed by our team.This innovative technique allows us to explore trapping phenomena within disordered systems confined to a finite volume.The three-body Padé approximant technique, central to our methodology, enables a nuanced understanding of the interplay between multiple excitations in a confined space.When extended to two-component systems, this technique reveals new facets of excitation dynamics that were previously inaccessible.Our study focuses on a finite volume system comprising immobile donor and trap molecules connected to a spherical nanoparticle.Excitations within this setup can be transferred between donors and from donors to acceptors (traps), with the trapping process assumed to be irreversible.
Furthermore, our approach provides a framework for exploring the effects of different spatial distributions and interactions among donors and acceptors.By varying parameters within our model, we can simulate different environmental conditions and predict how these changes impact the efficiency of energy transfer.This predictive capability is particularly valuable for designing experiments and interpreting experimental data in real-world applications.
Our methodology demonstrates remarkable versatility and can be applied to various molecular interactions.In this analysis, we specifically focus on dipole-dipole interactions.Our main objective is to determine the system's Green function, which allows us to compute all trapping and transport properties within the system.
The incorporation of multistep energy migration among groups of donors, followed by non-radiative excitation energy transfer to acceptors, is crucial for experimental measurements.These measurements require a sufficiently large number of initially excited donors to ensure accurate measurement of the fluorescent signal, considering potential disturbances and noise.The effectiveness of our model will be evaluated through Monte Carlo simulations.
Our Monte Carlo simulations not only corroborate the theoretical predictions but also highlight the practical applicability of our model.By demonstrating excellent agreement with the observed donor fluorescence intensity decay, we establish the credibility of our approach and its potential as a reliable tool for studying excitation transfer in various systems.

Theoretical model
Let us consider a spherical nanoparticle with radius R and surface S = 4πR 2 , containing N donor molecules denoted as D and M excitation energy traps denoted as A (acceptors), randomly distributed on the surface (refer to Fig. 1).The donor molecules are numbered from 1 to N , and the trap molecules are numbered from N + 1 to N + M .The individual configurations of molecules are described by the vector R = (r 1 , r 2 , . . ., r N+M ).
Let us introduce the following notation for the transfer rate for excitation energy migration w DD x i x j and w AA x i x j between two donors (D) and two acceptors (A) , respectively, while w DA x i x j represents the energy transfer rate from a donor to an acceptor.Moreover, self-transfer does not occur, so w DD x i x i ≡ 0 and w AA x i x i ≡ 0 .Each excited donor molecule, D * , may undergo deactivation through fluorescence emission (with rate constant k D F ), non-radiative transition (with rate constant k D q ) or non-radiative energy transfer from D * to another donor molecule D (with rate w DD x i x j ).The transfer rates k D F and k D q are related to the characteristic time τ 0D , by the relationship: The dynamics of electronic excitation energy migration and transfer in the system of spherical nanoparticles are described by the following master equation where P x i x j (t) denotes the conditional probability density of finding the excitation at time t on molecule x i given that at the initial moment t = 0, molecule x j was excited (where x i denotes all angular and spatial coordinates of the i − th molecule).The initial condition is P x i x j (0) = δ ij , where δ ij is the Kronecker delta.
The master equation can be expressed in matrix form as where (1) www.nature.com/scientificreports/ The solution to this equation is given by Information about the system is derived from the probability density, averaged over the spatial distribution R of molecules, that the excitation of the molecule at time t occurred at the point with the coordinate r where the bracket . . .R signifies the ensemble average over the molecules' distribution R.
The solution to the master equation is obtained using the Green's function method.The searched observable P(r , t ) is then expressed in the following form Let us assume that the initial probability distribution [P(R, t = 0)] j=1,N depends solely on the spatial distribu- tion of the excitation pulse, implying that only donor molecules were excited at the initial moment.
Then, the Green's function expressed by the elements of the matrix W assumes the following form The terms of the double sum above can be categorized into three classes: diagonal terms where j = k , and non-diagonal terms where j = k for j ≤ N and j > N .Thus the Green function G(r, r ′ ,t) can be represented as follows where The above Green functions can be interpreted as conditional probability densities of finding excitation at the point with the coordinate r and the moment of time t, provided that for t = 0 the donor was excited at the point with the coordinate r ′ .In view of further considerations, it is convenient to perform the Fourier-Laplace transform of the Green function G(r, r ′ ,t ).The transformed Green's functions are as follows: where the matrix E has the following elements The Green's functions G DD (k, ǫ), G DA (k, ǫ) and G SD (ǫ) are not independent.They collectively describe the probability density function, which should be normalized to unity The decay of the donor fluorescence is obtained by inverting the Laplace transform of the function:

Expansion in a series with respect to fluorophores concentration
In the following, we develop an expression for G D (k, ǫ) , the Fourier-Laplace transform of G D (t), which describes the donor fluorescence decay (see Eq. ( 16)).First, we calculate the Green function G DA (k, ǫ, N, M) , expanding it in a series in powers of the molecular surface density.For a system containing N donors and M acceptors distributed randomly over a sphere with a finite surface S , this expansion is as follows The term � DA 1 (k, ǫ, S) contains information about two-particle interactions (i.e.non-radiative excitation energy transfer).Similarly, � DA 2 (k, ǫ, S) and � DA 3 (k, ǫ, S) contain information about three-particle interactions, � DA 4 (k, ǫ, S) and � DA 5 (k, ǫ, S) four-particle interactions, and so on.If a general expression for the n − th term is found, G DA (k, ǫ, N, M) could be calculated exactly.In practice, only the first three terms can usually be evaluated.For a random distribution of molecules on the finite surface of a spherical nanoparticle, we need to average the above expression over all possible molecular positions.The explicit form of functions � DA i (k, ǫ, S) is presented in the Supplementary Material.

Padé approximant for donor fluorescence decay
Padé approximants are commonly employed in statistical mechanics to approximate truncated power series expansions, such as Eq.(17).Regardless of the number of terms calculated in the expansion for G DA (k, ǫ, N, M), the result may behave poorly for large N and M or small ǫ .In contrast, Padé approximants typically exhibit the correct asymptotic behavior in these limits and can serve as excellent approximations to G DA (k, ǫ, N, M).
By introducing the following functions where C(r 12 , r 13 , ǫ) and D(r 12 , r 13 , ǫ) are explicitly presented in the Supplementary Material, we obtain the fol- lowing expression for the Fourier-Laplace transform of the donor fluorescence decay G D (ǫ) The Padé approximant to G D (ǫ) is given by where f 2 (ǫ) > 0, f 3 (ǫ) > 0, and g 3 (ǫ) < 0. Considering the sign of the g 3 (ǫ) function (it is negative), it can be easily noticed, that including this term in the 3-body Padé approximant results in a slower fluorescence decay compared to the 2-body Padé approximant (see Eq. ( 26)).Conversely, inclusion of the f 3 (ǫ) term associated with energy migration in the donor set, which results in energy transfer to the acceptor in the 3-body Padé approximant, accelerates the donor fluorescence decay.The final effect depends on the balance between both effects.

2-body approximation
By introducing the variables the function f 2 (ǫ), representing interactions of two molecules (non-radiative energy transfer between one donor and one acceptor) is given by where a = 1 4 ξ 2 . By using the approximation described in the Supplementary Material (see Eq.(S21)), we obtain Thus, the 2-body Padé approximant of the Fourier-Laplace transform of the donor fluorescence decay is

3-body approximation
The function f 3 (ǫ), describing the interactions of three molecules (2 donors and 1 acceptor), when time is meas- ured in units of τ 0D with excluded natural decay exp − t τ 0D , assumes the following form (see Supplementary Material) where the integral J 3 (α, ǫ) can be approximated by Similarly, the function g 3 (ǫ) also describing the interactions of three molecules (1 donor and 2 acceptors), can be calculated as The expressions from Eq. ( 25), ( 27) and ( 29) inserted into Eq.( 22) finally give us the studied 3-body Padé approximant for G D (ǫ).
Figures 2a-d illustrate the results of donor fluorescence intensity decay in the presence of various numbers of acceptors.To provide a representative case for comparison, we selected a system comprising a nanoparticle with a radius R = 50 nm, labeled with donors and acceptors characterized by R DD 0 = 4.6 nm and R DA 0 = 5 nm 9 .The results in Figs.2a-d show the decay for different numbers of donors: 10 donors in Fig. 2a, 100 donors in Fig. 2b, 300 donors in Fig. 2c and 500 donors in Fig. 2d.This systematic variation allows us to evaluate the influence of energy migration preceding energy transfer on the donor fluorescence decay.Geometric symbols in the figures represent the results obtained from Monte-Carlo simulations, while the scattered lines depict the theoretical predictions within the three-body Padé approximant.It is evident from all figures that in the absence of acceptors, the predicted donor fluorescence decay follows a single exponential character and can be described by the same curve regardless of the number of donors.A similar result can be observed for three-dimensional solutions.However, as the number of acceptors increases, the fluorescence intensity of donors decays significantly faster, exhibiting a non-exponential character due to effective excitation trapping.Comparing the curves obtained for the same number of acceptors but different number of donors reveals that the donor intensity decays faster for systems with stronger energy migration.
The agreement between the Monte-Carlo result and the theory within the three-body Padé approximant is excellent up to 200 acceptors for small and moderate number of donors (Fig. 2a, b).Even in the case of a large number of donors, corresponding to strong energy migration, the agreement remains very good up to 100 acceptors.These results are surprisingly good considering the nature of the three -body Padé approximant described earlier.However, for an extremely high number of acceptors, the Monte-Carlo results predict faster decays compared to the theoretical predictions, reflecting the limitation of the Padé approximant.
During energy migration and transfer, excitation traverses many sites on a nanoparticle, making it crucial to understand the spatial extent of this stochastic process.A useful metric for this is the relative mean squared displacement of the excitation energy, denoted as (< r 2 >) 1/2 /R DA 0 .Figure 3 illustrates the Monte-Carlo results of the relative mean squared excitation displacement versus the number of donors for three fixed numbers of www.nature.com/scientificreports/acceptors.It is evident that the range of excitation walk increases significantly with the number of molecules present on the nanoparticle, reaching values between 1 and 1.35 depending on the number of acceptors.These values are substantial, indicating that excitation may reach distant sites from the initially excited one.Hence, energy migration should not be regarded as a purely local process in terms of nanoparticle size.
Figure 3 illustrates that in the case of a single donor (Förster case), the curves start higher, corresponding to a greater number of acceptors.This occurs because the emission of photons by the initially excited donor outweighs the energy transfer to the acceptors.As the number of donors increases, the migration of energy within the set of donors becomes increasingly crucial in determining the relative mean squared excitation displacement value.When the number of donors significantly exceeds the number of acceptors, energy migration becomes dominant.In such scenario, the curves reverse, with the curve corresponding to the smallest number of acceptors being the highest.This effect directly results from increased energy trapping with a growing number of acceptors, assuming a fixed number of donors.
Figure 4 depicts the results of Monte Carlo simulations of relative mean squared excitation displacement as a function of the number of acceptors for four fixed numbers of donors.In the absence of acceptors ( M = 0 ), the curves start higher, corresponding to a greater number of donors.This is attributed to the increasing importance of energy migration within the set of donors when energy transfer to acceptors is absent, leading to a higher relative mean squared excitation displacement with more donors.As the number of acceptors increases (with a fixed number of donors), the relative mean squared excitation displacement decreases due to increasing energy trapping.When the number of acceptors significantly exceeds the number of donors, the influence of energy migration diminishes, and all four curves asymptotically approach the Förster case.

Conclusion
The presented model introduces a novel perspective on energy migration and trapping on a spherical nanoparticle with chemically linked donors and acceptors.This model reveals that the dynamical behavior of the system is not only exquisitely responsive to the number of donors, acceptors, and the radius of the nanoparticle, but also showcases how these factors interplay in novel ways.Monte Carlo simulations, conducted for this innovative model, not only affirm the utility of the three-body Padé approximant but also highlight its novelty in capturing complex interactions.
Intriguingly, certain deviations observed between theoretical predictions and results of the Monte Carlo simulations at long times originate from the inherent characteristics of the Padé approximant, shedding new light on its limitations.While remarkably effective for scenarios involving a relatively modest number of attached fluorophores and short post-excitation durations, the Padé approximant's performance diminishes over extended timeframes and higher fluorophore counts.Nonetheless, its viability, even under demanding circumstances, remains noteworthy.
Our presentation of a reasonable analytical approximation holds considerable novelty, enabling facile prediction of spherical nanoparticle properties amidst energy migration and transfer.The high accuracy consistently verified through Monte Carlo simulations in numerous cases emphasizes its innovative contribution to near-"ideal" solutions.This underscores the substantial potential of leveraging Monte Carlo simulations, not only for simple spherical systems but also for more intricate core-shell nanoparticle configurations marked by higher complexities, where conventional analytical models often yield results that are only approximations.
In particular, our results demonstrate that the Monte Carlo technique can be effectively extended to model three-and four-component fluorescent systems on a nanosphere.This opens up possibilities for achieving a fully tunable system over the visible spectrum.The ability to accurately simulate and predict the behavior of such complex systems holds significant promise for the design and optimization of advanced photonic and optoelectronic devices.
Future studies will delve deeper into these multi-component systems, aiming to refine our understanding of energy transfer dynamics in even more intricate configurations.By continuing to leverage the Monte Carlo approach alongside the three-body Padé approximant, we anticipate uncovering further insights into the interplay of multiple excitations and their interactions in confined spaces.These efforts will pave the way for the development of new materials and technologies that exploit these sophisticated energy transfer mechanisms.
In conclusion, this study makes a significant contribution to the theoretical understanding of excitation transfer and trapping in disordered systems.The novel application of the three-body Padé approximant technique to two-component systems, validated by Monte Carlo simulations, opens up new avenues for research and development in materials science and related fields.Future work will focus on further refining the model, exploring its applications to other types of disordered systems, and extending the theoretical framework to include additional complexities such as varying donor-acceptor interactions and more intricate spatial configurations.

Figure 1 .
Figure 1.Geometrical relations between donors (D) and acceptors (A) on a spherical nanoparticle.

Figure 3 .
Figure 3. Monte Carlo results of relative mean squared excitation displacement as a function of the number of donors for three fixed numbers of acceptors.

Figure 4 .
Figure 4. Monte Carlo results of relative mean squared excitation displacement as a function of the number of acceptors for four fixed amounts of donors: N = 10;100;300;500.